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Abstract 

J> . We discuss the hadroproduction of charmed mesons in the framework of the 

constituent cascade model taking into account the valence quark annihilation. It 

$_i ' is shown that the small valence quark annihilation process dominates the leading 

particle production at large Feynman x and explains the recent experimental data 
on the asymmetry between D° and D° at 350 GeV/c. 
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1 Introduction 



Recently experiments at CERN [|TJ measured the neutral D mesons in 7r~ nucleus 
collisions and observed much smaller values of the leading/non- leading asymmetry than 
those of charged D mesons i.e. less than 0.2 and even a negative value around x = 0.8. The 
experimental data on the asymmetry of charged D mesons increases from zero to nearly 
one with Feynman variable x in the ir~ fragmentation region [p], 0, |3], |3|, || . The leading 
particle contains the same type of quark as one of the valence quarks in the incident 
hadron, while the non-leading one does not contain the projectile valence quarks. For 
example, the asymmetry of D° (cu) / D° (uc) in n~(du) interaction with nucleon is defined 
as 



A n - N (D°,D C 



a(D°) - a(D°) 
cr(D°) + cr(D°) 



In the perturbative QCD at leading order, the factorization theorem predicts that c 
and c quarks are produced with the same distributions and then fragment independently. 
In this case the asymmetry A n -^(D~ , D + ) is equal to zero ||. Even in the case of next to 
leading order, the predicted asymmetry is much smaller than the data JF]. The asymmetry 
A n ~^(D ~, D + ) has been investigated and explained by means of many approaches: string 
fragmentation [| 
process [15], [16|, 
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intrinsic charm contributions [§, [Hj, 
recombination using valon concept 



11, 12, 13, 14 



[0,0 



recombination 
and so on. We have 

proposed the constituent quark-diquark cascade model and explained the leading/non- 
leading asymmetry of charged D mesons A n - N (D~, D + ) successfully [fZIJ. The model, 
however, gives rather large values of A^-^^D , D°) at 0.5 ~ x as expected from the leading 
particle effect but deviating from the experimental data. Although several models have 

p(D°, D°) at x ~ 0.6 [III 17], the asymmetry 



been applied more or less satisfactorily to A^-p^D^, D°) at x ~ 0.6 
problem about the charmed hadron productions is an open question. 

In the present paper we investigate the leading/non- leading D-meson asymmetry in 
the framework of the constituent quark-diquark cascade model by taking into account the 
valence quark annihilation. 



2 Model description 

We consider an inclusive reaction A + B ^ C + X in the centre of mass system of A and 
B. The light-like variables of A and B are defined as follows: 

a E A ± p cm B _ E B p cm 



where ^J~s§ is the centre of mass energy of the incident hadrons A and B. We briefly 
review our model and then introduce the valence quark annihilation mechanism into the 
model. 
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Figure 1: The interaction mechanism in AB collision: a Non-diffractive dissociation type, 
b, c Single-diffractive and d double diffractive dissociation type mechanisms. 

When the collision between A and B occurs, the incident hadrons break up into two 
constituents with a probability (1 — P g i); otherwise they emit wee gluons with P g i followed 
by a quark-antiquark pair creation. We assume four interaction types: a) non-diffractive 
dissociation, b) and c) single-diffractive dissociations of B and A, and d) double-diffractive 
dissociation types as shown in Fig. |l|. The probabilities of these types to occur are 
(1 — P g i) 2 ,P g i(l — Pgi),P g i(l — P g [) and P^, respectively. Here we denote the quark- 
antiquark pair emitted from A (B) via the wee gluons as Ma ( Mb). The probabilities of 
M4 ( Mb) to be uu, dd, ss and cc are denoted as P ua , P^j Pss and P c5 , respectively. 

The momentum fraction of Ma is fixed by the distribution function 

H Ma/a (z) = z^-\\ - zf^ l /B{(3 gl ,(3 ld ), (3) 
and the uniform distribution R in the interval from zero to one as, 

Then the incident particles A and B have the following momentum fractions: 

x+ = Xq + (1-z), x a = m 2 A /(x^ s ), 



x_ — Xq_ (x_ Xq_(1 -R)), x + — Xq + , (5) 

where the mass shell condition is considered and transverse momenta are neglected. The 
momentum fraction of Mb is treated similarly, exchanging the role of A and B. 
In the centre of mass system of incidents A(Ma) and B(Mb), we define the light-like 
fractions of these hadrons and fix the light-like fractions of the projectile constituents. 
The distribution functions of the constituents in the projectile A composed of a and a' 
are described as 

Ha/A(z) - H a , /A (l - 2) = ±^-± . ,6) 
Then the light-like fractions of a and 

x a _ = Xq_R, x^_ = Xq + — x\ and x a _ = Xq_ — x a _, respectively. The distribution functions 
of the constituents in Ma, B and Mb are similarly defined. 



In our model hadrons are produced on the chain between a valence quark(anti-quark) 
from A(Ma) and the valence diquark(quark) from B(Mb) via the cascade processes 

q -> M(qq') + q', 



B(q[q'q"}) + [q'q"], B(q{q'q"}) + {q'q"}, 



[q'q"} -> B(q[q'q"]) + q, 



M(qq') + [qq"],M(qq') + {qq"}, 



{q'q"} -> £?(<?{<?'<?"}) +g, 



M(gg , ) + [?g"] 1 M( 9g -') + {g g "}, 

(7) 

where [ ] and {} denote the flavour antisymmetric and symmetric diquarks, respectively 



|2lfl . Meson production probabilities from q, [q'q"] and {q'q"} are 1 — e, ] and ?7{}, re- 
spectively. 

We redefine the light-like fractions of the incident constituents in the rest frame of the 
cascade chain. The momentum sharing of the cascade process q + q' — > M(qq") + q" + (f 
from a g with x± and q' with a;± takes place as follows p2 , 23 : First, using the emission 
function 

=z#<-\l-z)* + l , <t'- 1 /B(<yl3 q , (3 q + /V), (8) 

we fix the lightlike fractions of q" and Masi' = x+z and x¥ = x% — x% , respectively 
and put x q _ = x q _. Second, the transverse momentum of M is determined from the 
probability function 



Gd4) = ^«p(--^i4) (9) 



in space. Then, from the onshell condition, x^ is fixed as x^ 1 = {m 2 M + p T 2 ) / x+ s' , 
where y/s' is the subenergy of the cascade chain. The transverse momentum of q" is 
Pt = Pt ~~ Pt- The lightlike fraction of q' is decreased to xL = x q _ — x^ 4 . If the energy of 
q' is enough to create another hadron, the cascade such as q" + q' — > g" + g'" + M(q"'q') 
takes place in the opposite side. Finally recombined hadrons are put on-shell by a two 
body decay process as explained in [p0[|. 

The dynamical parameters /3's in (|6|) and (^), which determine the momentum sharings 
of the constituents, are related to the intercepts of the Regge intercepts as (3 U = Pd = 
1 - aip-ufi), (3 S = 1 - aty(0),/3 c = 1 - aj/^,(0) dH, |5f. From previous analyses |2|, |||, 
we determine the values for diquarks as /^-j = 7[ ]($ + /3j), /3^jy = 7{}(A + /3j)- We 
consider lower lying hadrons: pseudoscalar {PS), vector (V), tensor (T) mesons, octet (O) 
and decuplet(-D) baryons composed of u, d, and s flavours and the correspondings with 
charm flavour. We assume the production probabilities for them to be Pps, Pv-i Pt (= 
1 — P PS — P V ),P and Pjj (= 1 — Po), respectively. Octet and decuplet baryons are 
described as 

|8 >= cosO\q[q'q"] > +sm6\q{q'q"} >, (10) 
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Figure 2: The valence quark annihilation processes a soft annihilation and b semi-hard 
annihilation processes in ir~p collision. 

|io>=|<?W}>- (ii) 

Directly produced resonances decay into stable particles. Details of our model are ex- 
plained in 0, [ZT[ gg. 

We modify our model to include valence quark annihilation in the non-diffractive 
dissociation type interaction. Let us consider the case of n~p collision. We take account 
of the annihilation of the valence u from 7r~ and valence u from proton target in a slightly 
different way from the one pointed out in [|15|, 17]. There are annihilation processes 



such as shown in Fig. 2. The annihilation process in Fig. 2a is considered soft process 
and its contribution is related to the magnitude of a? nel — af nel by unitarity. This process, 
however, may be negligible for the charm-pair production. Here we only take into account 
the annihilation process in Fig. 2b and assume that the process 

uu — > qq (12) 

occurs with the probability P ann i(l — P g e) 2 and the non-diffractive type occurs with the 
probability (1 — P a nni)i)- — Pgt) 2 - Branching ratios of uu — > uu, dd, ss, and cc are chosen 
to be equal to each other for the channels allowed energetically. 

This process is considered as a semi-hard scattering process and produced q and q are 
supposed to be non-free. It seems natural to assume that q has tendency to be produced 
in the forward direction of the u in the centre of mass system of u and u, due to the 
confinement force between the valence d quark and the produced q quark. Here we choose 
the distributions for z = cos <fi as 

D(z) = 3 -(1 + Z y (13) 

in the region — 1 < z < 1, where <fi is the angle between the directions of u and q in the 
centre of mass system of u and u. After this annihilation mechanism, the non-diffractive 
type production of hadrons occurs as shown in Fig. 2b. 



3 Comparison with the data 

In this section, we give the results of our model for the inclusive hadron productions 
in 7i~N collisions. The value of the dynamical parameter (3 C is changed from the value 
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P c = 8.0 used in |20| to f3 c = 4.0. This is in accord with the argument that the slopes of 
Regge trajectories of charmed mesons are smaller than those of light mesons as discussed 
26 1 . The parameter of p T distribution C = 1.2 and the probability P c5 = 0.00016(P s j = 



in 



0.09984) in are changed into C = 1.0 and P c5 = 0.00025(P SS - = 0.09975). Furthermore 
the meson production probabilities are assumed as Pps = 0.4, Py = 0.4 and Pt = 0.2 for 
light mesons and Pps = 1/9, Py = 1/3 and Pt = 5/9 for charmed mesons. We put the 
parameter P an ni so as to reproduce the asymmetry A n - N (D°, D°). We choose the value 
Panni = 0.0005. This value scarcely changes the features of the spectra of light hadrons 
with u, d and s. For other parameters, we use the same values used in the previous analysis 

In Fig. 3 we show the results of D* productions and compare the two cases: case(l) 
3(1 + z) 2 /8 distribution for cos a and case(2) without annihilation (P an ni = 0:the same as 
in @ except for the values of (3 C = 4.0, C = 1.0, P cg = 0.00025, and P PS = 1/9, Py = 1/3 
and Pt = 5/9 for charmed meson productions). The annihilation effect of Fig. 2b is seen 
in D*°(uc) and D*~(dc). Our model gives a satisfactory description of x dependence for 
production of D* + or D*~ in tt~N collision. The result for production of D*° or D*° at 
x ~ 0.1 is small as compared with experimental data [27 |. 

The annihilation effect on meson productions in n~p is twice as much as in Tc~n colli- 
sion. The annihilation mechanism have a considerable effect on the non-leading particle 
D*°{uc). However, the shape of the leading particle D*~(dc) is affected little by the anni- 
hilation mechanism. Therefore our model gives larger differences between A n - p (D*°, D*°) 
and A w - n (P>*°,P*°) than those between A n _- p (D*- , D* + ) and [A n - n (D*~, £>*+). Fig.4a 
shows the results in case(l) for A n - P (D*°, D*°), A w - n (D*°, D*°) and A^-n^D* , D*°). 
Hereafter we show the average results of 7r~ beam on proton and neutron targets. Fig. 4b 
shows the results of A n - N (D*°, D*°) and A n - N (D*~, D* + ) in case(l) and case(2). 

Smaller values of A W - N (D*~, D* + ) as compared with A 1T - N (D* a , D*°) at x ps 0.5 are 
due to the difference between the hadron productions on quark-diquark and antiquark- 
quark chains in non- diffract ive dissociation type mechanism (Fig. la) ||20|| . The leading 
particle D*~(dc) is produced on the chain between the valence d quark in the beam and 
the valence diquark in the target. The proton target has a tendency to break into an 
energetic valence diquark and a wee valence quark. Then the D*~ meson produced in 
the first cascade step is energetic. However, in the case in which a baryon is produced 
from the valence diquark in the first cascade step, the momentum of the valence d quark 
is decreased. Furthermore the total momentum is shifted to the diquark side. Then the 
leading spectrum of D*~ on the quark-diquark chain tends to have a small momentum 
and the asymmetry is reduced in 0.3 ~ x ~ 0.7. 

On the other hand, the leading particle D*°(cu) is produced on the chain between the 
valence u quark in the tt~ beam and the valence quark in the target. The momentum 
reduction of the u quark is small even in a case a meson is produced in the first cas- 
cade step from the valence quark in the target fragmentation. Therefore the asymmetry 
A 1T - N (D*°, D*°) is not so suppressed as compared with A n - N (D*~, D* + ) in 0.3 ~ x ~ 0.7. 
In case(l), the annihilation effect compensates the quark-diquark effect and the asymme- 
try A^-n^D* , D*°) becomes negative at large x. 
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In Fig. 5, we show the result of x dependence of D mesons in ir~N collisions for case(l) 
and case(2). The annihilation effect is also seen in D° and a little in D~ spectra as in D* 
spectra. In Fig. 6, we show the result of p\ dependence of D mesons in tt~N collisions for 
case(l) and case(2). Our calculation shows smaller values of D° than those of D° at large 
p\ deviating from the data. The newly introduced annihilation process explains the large 
p\ charmed meson productions in part. Our model is in good agreement with the data 
except for D°, D° at large p\. In Fig. 7 the results of x dependence of A n -^(D^, D + ) 
and A n - N (D°, D°) are compared with the experimental data at 350 GeV/c |TJ. The 
agreement of x dependence with the experimental data is satisfactory. The negative value 
of A v - n(D°, D°) at large x is well explained by introducing a small amount of the valence 
quark annihilation process to the model. Fig.8 presents the results of p\ dependence of 
A 7T - N (D-,D + ) and A„- N {D Q ,D°). The p\ dependence of A n ~ N (D°, D°) disagrees with 
the data in sign as seen in Fig.6b. The discrepancy increases with the annihilation effect. 
The results of the x and p\ dependences of A n ~ N (D~, D + ) and A n - N (D°, D°) at 500 
GeV/c are compared with the experimental data @| in Fig.9. The agreement with the 
data is fairly good except for the negative value of the data at x ~ —0.15. There is 
a valence d quark in the incident nucleon and our model predicts a positive value of 
A n - n{D~ , D + ) for x < as seen in Fig.9. However, in the case of proton target with 
cos 6 = 1 in (|10"1), there is no valence d quark and the model gives a small value of 
A^-xiD^ , D + ) in the target fragmentation region. 

The results of the x and p\ dependence for production of Dg(sc) or Dg(cs) are 
compared with the experimental data |I| in Fig. 10. The features of the spectra are in 
good agreement with the data except for the small discrepancy in normalization. Fig. 11 
shows the calculated results of asymmetries A w - N (Dg, Dg) with respect to x and p\. 
The values are small because both mesons are non-leading particles. But in case(l) 
A n -j^(Dg , Dg) increases with x at large x due to the valence quark annihilation. 

4 Discussions 

The large leading/non- leading asymmetry A n -N(D~, D + ) is naturally explained by our 
constituent quark cascade model as in the case of light hadron productions. In [|(J , we no- 
ticed that the cascade chain properties are different between the antiquark-quark and the 
quark-diquark chains. The leading particle D~(dc) produced in the quark-diquark chain 
in 7i~p collision is less energetic when a baryon is produced in the target fragmentation in 
the first cascade process as compared with D + (cd) produced in the quark-antiquark chain 
in 7i + p collision. We have different leading particle effects between 7r + and 7r~ beams i.e. 
A W + N (D + , D~) > A n - N (D~ , D + ) around x ~ 0.6 PD| . Although the difference between 
the asymmetries A T - N (D°, D°) and A T - N (D~, D + ) disappears due to the decay effect 
and the valence quark annihilation effect, this quark-diquark chain effect is seen in the 
difference between the asymmetries A n -(D*°,D*°) and A n - N (D*~, D* + ) 
around x ~ 0.5 as shown in Fig. 4b. 

The small valence quark annihilation process (^) explains the negative value of 
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A W - N (D°, D°) at large x observed in the recent experiment |T|]. This process is considered 
as semi-hard interaction. All channels occur with equal probability and the produced 
constituents tend to have the forward direction of the incident valence constituents. This 
is a part of the hadron production process at large p\. Our model can explain the spectra 
up to the region p\ ~ 10 (GeV/c) 2 . As noticed in the previous section, the results of cal- 
culated behaviours of the p\ dependence of ^-^(-D , D°) disagree with the experimental 
data. This implies that there may be considerable contribution from soft interactions 
which maintain the leading particle effect even at large p\ region. The negative values 
of experimental data on A^-nID , D°) at p\ ~ 3 (GeV/c) 2 suggests that the q quark in 
the annihilation process ( |T2D is produced in the very small cone around the direction of 
incident u as compared with flllf ). It is interesting to investigate the leading/nonleading 
asymmetry in KN collisions, since in K + N collision, there is no annihilation process. 
These points will be discussed elsewhere. 
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Figure 3: Comparison of the model with the experimental data on the x-dependences for 
productions of a D* + or D*~ and b D*° or D*° mesons at P L = 360 GeV/c J2?|. The 
theoretical lines were calculated for Pi = 350 GeV/c, full line for case(l) with the valence 
quark annihilation and dotted for case(2) without the valence quark annihilation. 

Figure 4: a The x- dependences of the asymmetries A n - p (D*°, D*°), A w - n (D*°, D*°) and 
A^-n^D* , D*°). b The comparison of x-dependence of the asymmetries A^-^iD* , D*°) 
and A T - N (D*~ , D* + ) between in the case(l) and case(2). 

Figure 5: The x-dependences of a D and b D° and D° mesons at Pi = 350 GeV/c. The 
full lines were calculated for case(l) with the valence quark annihilation and dotted for 
case(2) without the valence quark annihilation. The experimental data are taken from 

Figure 6: Same as in Fig.5 for p^-dependence in range < x < 1. 



Figure 7: The x-dependences of asymmetries a A^-n^D^, D + ) and b A^^{D , D°) at 
Pi = 350 GeV/c. The full lines show the results in case(l) and dotted lines in case(2). 
The experimental data are taken from 

Figure 8: Same as in Fig. 7 for p^-dependences in range < x < 1. 

Figure 9: The asymmetry A n - N (D~, D + ) at p L = 500 GeV/c a as a function of x, b as 
a function of p\ for range of —0.2 < x < 0.8 and 0.4 < x < 0.8. The experimental data 
are taken from ||]. 

Figure 10: Differential cross sections with respect to a x and b p\ for production of D\ 
or Dg at Pi = 350 GeV/c. The full line denotes the results for case(l) and dotted for 
case (2). The experimental data are taken from |]J. 

Figure 11: Asymmetries A n - N (Dg , D$) with respect to a x and b p\ for at Pi = 350 
GeV/c. The full line denotes the results for case(l) and dotted for case(2). 
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